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ABSTRACT 


Preliminary  results  of  an  experimental  investigation 
of  heat  transfer  to  liquid  helium  under  forced  flow 
conditions  are  reported  for  a  0.  213  cm  i.  d.  x  10  cm 
long  test  section  subject  to  the  following  range  of 
operating  conditions: 

System  pressures  1.1  -  2.1  atm  ^ 

Mass  velocities  4-64  g/s-cm 

Heat  fluxes  0 .  04  -  0 .  53  W/cm^ 

Inlet  subcooling  0.03  -  0.  10  K 

The  effect  of  the  above  system  parameters  on  the  heat 
transfer  and  critical  heat  flux  is  discussed;  a  compari- 
son of  forced  convection  boiling  with  other  modes  of 
heat  transfer  (pool  boiling  and  supercritical)  and  the 
performance  of  a  centrifugal  pump  used  for  circulating 
the  liquid  helium  are  also  included  in  the  report. 

Key  Words:    centrifugal  pump;  critical  heat  flux; 
film  boiling;  forced  convection;  heat  transfer;  helium; 
nucleate  boiling;  subcritical;  supercritical. 
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FORCED  CONVECTION  HEAT  TRANSFER  TO 
SUBCRITICAL  HELIUM  I 


1.  INTRODUCTION 

In  a  previous  study  by  Giarratano,  Arp  and  Smith  [l]  at  the 
Cryogenics  Division  of  NBS,  Boulder,  Colorado,  heat  transfer 
coefficients  for  forced  flow  of  supercritical  helium  were  measured 
and  a  correlation  was  developed  to  predict  the  heat  transfer  in  this 
region.    As  an  extension  of  that  study  this  report  describes  the 
preliminary  experimental  results  obtained  for  forced  convection 
heat  transfer  in  the  subcritical  region,  i.  e.  ,  below  the  critical 
pressure  2.  245  atm.*  Operation  of  the  centrifugal  pump,  used  for 
circulation  of  helium  around  the  test  loop  in  the  subcritical  region 
will  also  be  discussed. 

2.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 
AND  MEASUREMENT 

A  schematic  of  the  boiling  heat  transfer  flow  loop  is  shown  in 
figure  1.  A  centrifugal  pump,  previously  described  in  [2],  circu 
lates  the  liquid  around  the  flow  loop.  To  maintain  a  constant  inlet 
temperature  to  the  test  section,  heat  is  removed  from  the  liquid  a 
it  passes  through  approximately  13  meters  of  0.  32  cm  i.  d.  copper 
tubing  located  in  a  heat  exchanger  reservoir  of  approximately  2.  5 
liter  capacity.  The  heat  exchanger  reservoir  is  continuously  sup- 
plied with  liquid  helium  from  a  storage  dewar. 


•j"  Identification  of  any  materials  or  their  manufacturer  by  the  National 
Bureau  of  Standards  in  no  way  implies  a  recommendation  or  endorse- 
ment by  the  Bureau.    Furthermore,  use  of  other  trade  names  in  this 
paper  is  for  the  sake  of  clarity  and  does  not  in  any  way  imply  a  recom- 
mendation or  endorsement  by  the  Bureau. 
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*  1  atmosphere  =  1.  013  x  10    MN/m  . 


2.  1  Test  Section 

The  test  section  is  a  0.  213  cm  i.  d.  x  20  cm  long  stainless 
steel  tube  with  a  wall  thickness  of  0.  016  cm.    It  is  resistance 
heated  along  10  cm  of  its  length  (L/D  =50).    A  stainless  steel- 
pyrex  seal  located  in  the  bottom  of  the  loop  provides  electrical 
isolation  of  the  test  section.    A  wire  wound  preheater  upstream 
of  the  test  section  allows  variation  of  quality  at  the  inlet  of  the 
tube.    Accurate  measurement  of  the  current  and  potential  drop 
across  the  two  heaters  is  used  in  the  calculation  of  heat  flux. 
At  the  lowest  power  levels  the  uncertainty  is  less  than  2%. 

2.  2  Temperature  Measurement 

Outside  tube  wall  temperatures  are  measured  at  ten  points 

along  the  heated  test  section  length  with  calibrated  carbon 

resistance  thermometers.    A  detailed  sketch  of  a  thermometer 

t 

and  the  mounting  scheme  is  shown  in  figure  2.      The  thermome- 
ters are  positioned  one  cm  apart  with  the  first  thermometer 
located  0.  6  cm  from  the  inlet.    Inside  wall  temperatures  are 
obtained  by  calculation  of  the  temperature  drop  through  the  tube 
wall  using  the  thermal  conductivity  of  stainless  steel  given  by 
equation  (Al)  in  the  appendix  of  reference  [l].    Bulk  fluid  tem- 
peratures are  measured  at  four  points  along  the  flow  loop 
(upstream  of  preheater,  upstream  of  test  section  and  two  meas- 
urements downstream  of  the  test  section)  with  germanium 
resistance  thermometers.     These  thermometers  are  plotted  with 
vacuum  grease  in  copper  wells  which  are  soldered  to  the  tube. 


2 


3 


Electrical  Leads  Tempered 
on  Copper  Clamp 

Copper  Clamp 


Radiation  Shield  Thermally 
Anchored  to  Clamp 

Carbon  Resistance  Thermometer 
with  one  Lead  Soldered  to  Clamp 

0.0127  cm  thick  (.005")  Mylar 


Test  Section 


Figure  2. 


'Sketch  of  mounting  scheme  for  carbon  resistance  thermometers. 


The  carbon  resistance  thermometers  and  the  germanium 
resistance  thermometers  used  in  this  study  were  calibrated  in 
a  separate  apparatus,  over  the  range  4  -  20  K,  against  three 
germanium  resistance  thermometers  (GT  1011,  1027  and  1024), 
These  three  germanium  thermometers  had  been  previously- 
calibrated  against  NBS  secondary  standard  germanium  ther- 
mometers (GT  722  and  734). 

In  our  calibration,  over  the  range  4  -  9  K  the  agreement 
between  GT  1011,  1027  and  1024  was  within  0.  OIK  for  the 
worst  case  and  was  generally  within  0.  005  K.    Above  9  K  the 
agreement  was  within  0.  05  K.    An  average  value  of  the  three 
germanium  thermometer  readings  was  taken  to  be  the  true 
temperature. 

The  calibration  data  for  each  thermometer  were  fitted  with 
a  curve  of  the  form 

LOG  T  =  r  A(N)  LOG 

with  an  rms  percent  error  of  0.  300  in  temperature  or  less. 

After  a  period  of  about  six  months,  and  after  several 
intermittent  thermal  cycles  from  room  temperature  to  helium 
temperature,  three  carbon  thermometers  and  two  germanium 
thermometers  were  removed  from  the  experimental  apparatus 
and  dip  tested  in  liquid  helium  (4.  025  K)  to  check  for  a  shift 
in  calibration.    The  carbon  resistors  were  within  0.  1%,   1.  1%, 
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and  1.  2%  in  temperature  of  their  original  calibration  while  the 
germanium  thermometers  were  within  0.  3%  and  0.  5%,  It 
should  be  noted  that  while  this  does  not  represent  a  calibration 
shift  for  the  germanium  thermometers^  the  error  for  two  of  the 
carbon  resistors  lies  slightly  outside  of  the  99%  confidence 
limit  of  the  original  calibration  (1,  e.   3  X  rms  percent  error). 

In  the  two-phase  region,  since  the  pressure  drop  across 
the  test. section  was  of  the  order  of  a  few  mm  Hg  for  the  data 
presented,  the  intermediate  bulk  temperatures  were  taken  to 
be  the  saturation  temperatures  corresponding  to  the  inlet  static 
pressures.    This  approximation  was  within  experimental  error. 

Periodic  zero  power  runs  were  made  to  check  consistency 
among  the  fourteen  thermometers.    Agreement  between  the 
four  germanium  resistance  thermometers  was  within  a  few 
millidegrees  while  the  carbon  resistance  thermometers  showed 
a  decreasing  offset  with  time  when  compared  to  the  germanium 
resistance  thermometer  readings.    The  time  constant  (of  the 
order  of  hours)  was  too  long  for  it  to  be  practicable  to  wait. 
Since  the  time  of  each  powered  run  and  each  zero  power  run 
was  recordedjthis  well-defined  offset  was  used  to  correct  the 
carbon  resistance  thermometer  readings.    This  correction 
alsb  takes  care  of  any  shift  in  calibration  of  the  carbon  resist- 
ance thermometers.    E}xtensive  tests  on  the  carbon  resistance 
thermometers  in  a  separate  apparatus  (with  the  thermometers 
mounted  exactly  the  same  as  in  the  present  study)  showed  that 
they  are  extremely  sensitive  to  irradiation  from  high  tempera- 
ture sources.     It  is  felt  that  although  substantial  effort  was 
devoted  to  radiation  shielding  of  the  carbon  thermometers,  the 


6 


radiation  could  not  be  totally  eliminated  during  the  period  that 
the  entire  apparatus  approached  temperature  equilibrium. 
The  estimated  uncertainty  in  outside  wall  temperatures  is  at 
most  0,  05  K.    The  error  in  the  temperature  difference  between 
the  inside  wall  and  the  bulk  temperature  is  due  to  uncertainty 
in  the  thermal  conductivity  of  the  stainless  steel  wall,  the  wall 
thickness,  and  uncertainty  in  bulk  temperature.    This  varies 
from  26%  to  53%  below  the  critical  heat  flux  (see  3.  1)  to  0.  3% 
above  this  transition. 

2.  3  Pressure  and  Flow  Measurement 

Provision  was  made  to  measure  pressure  drop  across  the 
flow  orifice,  preheater,  and  test  section, and  the  static  test 
section  inlet  pressure  using  calibrated  pressure  transducers 
and  a  static  Bourdon  tube  pressure  gauge  accurate  to  0.01  atm 
located  outside  the  cryostat  at  room  temperature.    Test  sec- 
tion pressure  drop  will  be  reported  in  a  future  report  but  was 
small,  generally  less  than  7  mm  Hg. 

Since  the  fluid  at  the  outlet  of  the  preheater  was  slightly 
subcooled  for  the  data  presented  here, a  calorimetric  method 
determined  the  flowrate.    The  calorimetric  flowrate  was  also 
used  to  determine  a  discharge  coefficient  for  the  flow  orifice 
so  that  the  orifice  may  be  used  for  flow  measurement  when 
calorimetric  flow  determination  is  not  possible,  i.  e.  ,  two- 
phase  out  of  the  preheater.    The  discharge  coefficient  so 
obtained  agrees  well  with  that  of  an  identical  orifice  section 
which  was  calibrated  in  a  separate  apparatus.    The  accuracy  of 
the  flow  rate  measurement  has  not  yet  been  fully  evaluated,  but 
we  believe  the  values  based  on  calorimetry  to  be  accurate  to  a 
few  percent. 


7 


2.4  Extraneous  Heat  Exchange 

To  minimize  heat  leak  from  room  temperature  to  the  test 
section  by  conduction,  all  electrical  leads  and  the  thin  walled, 
0.317  cm  diam.  stainless  steel  pressure  transmission  lines 
are  thermally  anchored  to  the  outside  copper  surface  of  the 
liquid  helium  heat  exchanger.    A  length  of  multifilament 
niobium-titanium  superconducting  wire  is  used  for  power  leads 
between  the  heat  exchanger  and  the  test  section  and  between 
the  heat  exchanger  and  the  preheater.    The  small  diameter  and 
low  thermal  conductivity  of  this  wire  further  minimizes  heat 
leak  due  to  conduction  and  joule  heating  in  the  leads  is  elimi- 
nated.   Error  in  heat  flux  due  to  axial  conduction  from  the  ends 
is  less  than  2%  for  the  worst  condition.    Extraneous  heat  to  and 
from  the  test  section  is  therefore  considered  negligible.  Boil- 
off  from  the  liquid  helium  heat  exchanger  is  routed  through 
coils  soldered  to  the  pump  housing  and  the  radiation  shield. 
This  arrangement  minimizes  heat  leak  from  room  temperature 
via  conduction  along  the  pump  housing  and  provides  a  low  tem- 
perature radiation  shield  around  the  flow  loop.    The  test  sec- 
tion portion  is  further  protected  from  radiation  by  a  copper 

shield  thermally  anchored  to  the  heat  exchanger.    The  evacu- 

-7 

ated  copper  enclosure  (vacuum  less  than  10      mm  Hg)  which  is 
submerged  in  a  bath  of  liquid  nitrogen  provides  first  stage 
radiation  shielding  from  room  temperature. 

2.  5  Experimental  Measurement 

For  a  fixed  system  pressure,  pump  speed  and  quality  at 
the  inlet  of  the  test  section,  the  power  to  the  test  section  was 
varied  over  a  range  from  zero  to  that  at  which  a  discontinuous 
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rise  in  the  wall  temperature  occurred  at  the  outlet  end  of  the  test 
section  (critical  heat  flux).    For  some  of  the  runs  this  heat 
flux  was  exceeded  causing  the  discontinuity  in  wall  temperature 
to  move  up  the  test  section.    After  a  change  in  test  section 
power,  thermometer  voltages  stabilized  generally  within  a  few 
seconds  and  were  recorded  by  a  digital  voltmeter  and  auto- 
matically punched  on  paper  tape  together  with  all  other  pertin- 
ent voltages  and  run  information. 

Since  the  flow  loop  is  a  closed  system,  as  power  was 
applied  it  was  necessary  to  vent  the  system  to  maintain  a  con- 
stant test  pressure, and  conversely  a  decrease  in  power  applied 
required  adding  and  condensing  gas  in  the  loop  to  maintain  the 
pressure. 

This  procedure  was  repeated  for  different  pump  speeds, 
inlet  quality  and  pressure. 

During  a  measurement  the  temperatures  and  pressure 
were  stable  to  within  their  accuracy  prior  to  a  transition  in  the 
heat  transfer  mechanism  (critical  heat  flux  exceeded).  When 
a  transition  occurred,  it  was  accompanied  by  fluctuations  in 
temperature  (for  the  wall  stations  in  the  transition  region)  of  a 
rather  random  frequency,  and  an  amplitude  of  the  order  of 
0.  5  K.    Oscillations  in  pressure  were  not  noticeable 
but  this  was  at  least  partly  accountable  for, the  pressure  trans- 
mission line  being  heavily  damped. 

One  of  the  goals  of  the  present  study  is  to  establish  the 
repeatability  of  a  measurement  in  the  form,  say,  of  the  spread 
in  wall  temperatures  which  can  be  expected  when  the  flowrate, 
pressure  and  inlet  enthalpy  are  specified.    To  date,  we  do  not 
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have  sufficient  information  to  establish  a  day-to-day  variability, 
but  we  do  have  a  few  figures  for  the  variability  observed  during 
the  course  of  a  day's  running.    Below  the  critical  heat  flux  this 
variability  is  generally  of  the  order  of  a  few  percent  of  the 
excess  of  wall  temperature  over  bulk  temperature,  with  a 
worst  value  observed  of  11%.    Above  the  transition,  we  do 
not  have  sufficient  data  at  this  time  to  estimate  the  variability. 

HEAT  TRANSFER  RESULTS 
3,  1  Wall  Temperature  Profiles 

Temperature  profiles  along  the  wall  for  subcritical  helium 

are  shown  in  figures  3  through  10  for  pressures  from  1.  1  to 

2 

2.  1  atm  and  mass  velocities  from  4.  8  to  63.  6  g/cm  .  A 
typical  profile  is  extremely  flat  up  to  some  point  at  which  a 
sharp  rise  is  observed  for  the  higher  heat  fluxes.    This  we 
identify,  with  reference  to  a  vast  body  of  literature  on  boiling 
heat  transfer,  as  indicating  a  hydrodynamic  transition  from 
wetted-wall,  e.g.  ,  nucleate  boiling,  to  dry-wall,  e.g.  ,  film 
boiling,  heat  transfer.    Thus,  at  a  particular  point  along  the 
test  section  there  is  a  critical  heat  flux  above  which  the  wall 
temperature  rises  quite  steeply  with  heat  flux. 

Below  the  critical  heat  flux,  heat  transfer  is  very  good  as 
evidenced  by  the  wall  temperatures  being  nowhere  more  than 
0.  3  K  above  the  bulk  fluid  temperature.    Indeed  this  behavior 
is  strongly  reminiscent  of  boiling  heat  transfer  without  forced 
convection,  i,  e,  ,  pool  boiling.    In  figures  11  and  12  we  plot 
the  heat  flux  q  against  the  difference  between  the  wall  and 
bulk  fluid  temperatures  for  the  first  and  last  wall  thermometers 
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Figure  3.    Wall  temperature  profile  for  test  sectioru 
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Pressure  =  1.1  atm 

Mass  Velocity  =  11.3  g/s  cm^ 

inside  Dia.  of  Test  Section  =  0.213  cm 


q  =  0.257  W/cml 
X  i  Outiet  =  0.23 


0.247 
0.21 
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q  -  0.180 
X  =  0.15 
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Figure  4.    Wall  temperature  profile  for  test  section. 


LENGTH  ALONG  TEST  SECTION,  cm 


Figure  5.    Wall  temperature  profile  for  test  section. 
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Pressure  =  1.2  atm 
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LENGTH  ALONG  TEST  SECTION,  cm 

Figure  6.    Wall  temperature  profile  for  test  section. 
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Pressure  =  1.34  atin 

Mass  Velocity  =  63.6  g/s-cm^ 

Inside  Dia.  of  Test  Section  =  0.213  cm 
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LENGTH  ALONG  TEST  SECTION,  cm 

Figure  7.    Wall  temperature  profile  for  test  section. 
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Pressure    1.5  atm 

Mass  Velocity   15.3  g/s  cm^ 

Inside  Oia.  ot  Test  Section   0.213  cm 


q  0.308  W/cm^ 
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Figure  10.    Wall  temperature  profile  for  test  section. 
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Figure  11,    Heat  flux  vs  temperature  difference  for  station  1. 
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Figure  12.    Heat  flux  vs  temperature  difference  for  station  10. 
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respectively.    These  are  at  3  and  48  tube  diameters  respec- 
tively from  the  commencement  of  the  heated  length  of  the  test 
section.    We  have  also  plotted  as  a  solid  line  the  nucleate  pool 
boiling  correlation  of  Kutateladze  [7]  and  indicated  at  its  upper 
extremity  the  critical  heat  flux  for  pool  boiling  given  by  a 
second  correlation  due  to  Kutateladze  [7  ]  .    As  shown  for 
example  by  Brentari  and  Smith  [8]  these  correlations  give  a 
good  average  representation  of  nucleate  pool  boiling  data  for 
cryogenic  fluids.    The  present  data  for  forced  convection  heat 
transfer  are  represented  by  the  first  of  these  correlations 
below  the  critical  heat  flux  as  well  as  any  given  set  of  pool 
boiling  data,  which  are  notorious  for  their  sensitivity  to  the 
precise  nature  and  preparation  of  the  boiling  surface,  as  well 
as  its  heating  history  (a  strong  hysteresis  effect  is  often 
observed).    We  conclude  that  forced  convection  has  had  very 
little  effect  on  the  rate  of  heat  transfer  at  heat  fluxes  below 
the  critical  and  that  the  boiling  mechanism  itself  is  the  pri- 
mary determinant. 

In  figure  13  a  wall  temperature  profile  for  subcritical 
helium  is  compared  with  one  for  supercritical  helium  flowing 
at  approximately  the  same  mass  velocity  and  the  same  applied 
heat  flux.    This  comparison  illustrates  well  the  essential  dif- 
ference between  two-phase  and  single-phase  heat  transfer; 
provided  transition  is  not  allowed  to  occur  the  boiling  phenome- 
non is  extremely  effective  in  holding  down  the  wall  temperature. 
However,  we  draw  attention  to  the  last  wall  temperature  which, 
in  the  subcritical  profile  shows  a  slight  increase  in  tempera- 
ture indicative  of  the  onset  of  transition  to  dry-wall  heat  trans- 
fer.   This  is  seen  in  figure  9  in  a  more  highly  developed  stage 
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Figure  13,    Comparison  of  subcritical  and  supercritical  wall  temperature  profiles 
under  similar  mass  velocity  and  heat  flux  conditions. 
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at  higher  heat  fluxes.    The  comparison  between  heat  transfer 
in  these  two  pressure  regions  is  discussed  further  in  3.  3 
below. 

.  2  Critical  Heat  Flux 

Beyond  the  critical  heat  flux  the  wall  temperatures  may- 
rise  to  quite  unacceptable  heights,  from  the  point  of  view  of 
the  cooling  of  superconductors,  as  is  clearly  seen  from 
figures  3  through  10  and  figure  12.    It  is  therefore  of  great 
importance  to  be  able  to  predict  the  conditions  under  which  the 
transition  from  wetted-wall  to  dry-wall  heat  transfer  takes 
place.    In  fact,  from  this  point  of  view  it  is  almost  the  only 
consideration  of  importance  until  the  flow  rates  become 
very  high,  because  the  wall  temperature  rises  observed  in  the 
wetted-wall  region  are  so  small  as  to  be  of  no  concern,  and 
in  the  dry-wall  region  they  are  generally  so  large  as  to  preclude 
this  region  from  consideration.    However,  at  high  flow  rates, 
as  seen  in  figure  7,  the  wall  temperature  excursion  appears 
to  have  reached  a  limit  which  may  be  of  some  practical  value 
in  superconducting  technology.    Figure  7  gives  data  for  the 
highest  flow  rates  obtained  in  this  study,  and  thus  it  appears 
that  at  very  high  flow  rates  the  transition  may  be  kept  within 
useful  bounds. 

One  of  the  ultimate  goals  of  this  program  will  be  to  provide 
a  means  of  predicting  the  transition,  that  is:    given  the  state  of 
the  fluid  at  entry  to  the  heated  section,  i.e.,  enthalpy  and 
pressure,  the  flow  rate  and  the  heat  flux,  we  wish  to  know  at 
what  position  downstream  the  transition  will  occur.  Alterna- 
tively, given  the  state  and  flow  rate,  we  wish  to  know  the  heat 
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flux--the  critical  heat  flux- -at  which  the  transition  will  occur 
at  a  certain  position.    Considerably  more  data  than  are  so  far 
available  from  this  study  will  be  required.    For  the  present, 
however,  we  shall  confine  our  attention  to  identifying  trends. 

The  first  and  most  obvious  trend  is  that,  for  a  given 
pressure  and  mass  velocity,  the  critical  heat  flux  is  a  function 
of  distance  from  the  inlet  to  the  heated  section.    In  experiments 
on  heat  transfer  to  subcritical  helium  under  conditions  of 
natural  convection,  Johannes  and  Mollard  found  that  most  of 
their  data  could  be  represented  by  a  single  curve  when  the 
critical  heat  flux  was  plotted  against  the  distance  from  the 
inlet  to  the  heated  section  in  diameters  (the  hydraulic  or  equiv- 
alent diameter   D     is  used  for  rectangular  cross  sections 
e 

where   D    =  4  x  cross-sectional  area/heated  perimeter), 
e 

Figure  14  shows  this  plot  with  our  data  points  superimposed. 
The  short  channel  data  of  James,  et  al,  ,  the  long  channel 
data  of  Keilin,  et  al,  ,  and  the  data  of  Jergel  and  Stevenson 
are  also  included.    Those  of  Jergel  and  Stevenson  were  for  a 
long  channel  with  a  small  heater  at  the  inlet.    These  and  the 
data  of  Keilin,  et  al.  ,  were  for  forced  flow. 

The  trend  shown  by  Johannes  and  Mollard  is  confirmed, 
but  now  also  an  enhancement  with  increasing  mass  velocity  is 
observed.    The  point  at  2.  1  atm  (reduced  pressure  0.  93)  being 
lower  even  than  the  curve  representing  the  data  of  Johannes 
and  Mollard  is  not  surprising  and  illustrates  a  fact  familiar 
from  pool  boiling,  namely,  that  the  critical  flux  goes  to  zero 
at  the  critical  pressure  in  accord  with  Kutateladze's  critical 
heat  flux  correlation. 
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Figure  14.    Critical  heat  flux  as  a  function  of  length  to  diameter  ratio  with  pressure 
and  mass  velocity  as  parameters. 
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It  is  hardly  necessary  to  point  out  that  if  the  conditions 
leading  to  transition  in  forced  convection  heat  transfer  were 
the  same  as  those  for  pool  boiling, then  the  entire  tube  would 
undergo  transition  at  the  same  heat  flux  contrary  to  figure  14. 
It  is  therefore  necessary,  in  order  to  properly  account  for 
transition  in  forced  convection  heat  transfer,  to  look  for  either 
purely  local  conditions  that  change  with   L/D   or  to  relate  the 
hydrodynamic  condition  at  a  given   L/D  to  the  entire  set  of 
conditions  upstream.    In  the  search  for  a  local  explanation  it 
is  tempting  to  plot  transition  data  against  the  local  mass  frac- 
tion of  vapor  (quality)  at  a  fixed   L/D.    Since  the  heat  transfer 
coefficient   (h  =  q/  AT)    is  a  sensitive  indicator  of  hydrodynamic 
conditions,  we  have  plotted  this  in  figure  15  against  quality 
at   L/D  =  48.    Transition  is  strongly  indicated  by  a  pronounced 
maximum.    We  note  the  following: 

i.  As  with  the  critical  heat  flux,  the  critical  heat  transfer 
coefficient  increases  with  mass  velocity  for  a  given 
pressure. 

ii.  The  transition  occurs  at  decreasing  qualities  for 
increasing  mass  velocity  for  a  given  pressure, 

iii.  For  a  given  quality  the  heat  transfer  coefficient  gener- 
ally increases  with  mass  velocity  below  the  transition, 

iv.  The  higher  the  quality  at  which  transition  occurs  the 
lower  the  corresponding  heat  transfer  coefficient. 

V.  For  the  pressures  plotted  in  figure  15,  increasing 
pressure  has  the  same  effect  as  increasing  mass 
velocity.    Unfortunately  we  do  not  have  sufficient  data 
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Figure  15.    Heat  transfer  coefficient  as  a  function  of  quality  with  mass  velocity  and 
pressure  as  parameters. 


27 


at  2,  1  atm  to  be  able  to  generalize  this  observation  to 
near-critical  pressures. 

vi.    In  all  cases  transition  occurs  while  a  substantial 
fraction  of  liquid  remains. 

It  is  clear  that  successful  correlation  of  transition  obser- 
vations presents  a  considerable  challenge;  there  are  at  least 
six  interdependent  variables  involved  apart  from  orientation 
effects.    These  are  the  critical  heat  flux,  mass  velocity,  quality, 
pressure,  diameter,  and  position  along  the  heated  tube  or  channel. 
We  know,  however,  that  this  correlation  must  reduce  to  that 
of  Kutateladze  at  low  mass  flow  rates,  for  then  the  transition 
can  be  understood  as  a  critical  ratio  between  inertial,  interfacial 
and  buoyancy  forces,  which  was  the  basis  of  that  correlation. 
The  development  of  a  suitable  correlation  will  be  left  to  a  sub- 
sequent phase  of  this  work. 

A  rather  interesting  comparison  can  be  made  for  the  pres- 
ent transition  results  and  those  of  Keilin,  et  al.  ,  for  similar 
mass  velocities  and  pressures.    These  authors  report  heat 
transfer  coefficients  in  flat  horizontal  coils  of  copper  tubing  of 
0.45  mm  inside  diameter  at   L,/D  =  278.    In  figure  16  the  com- 
parison is  made  with  our  data  at  similar  mass  velocities.  The 
significant  difference  is  in  the  quality  at  which  transition  appar- 
ently takes  place.    The  heat  transfer  coefficient  is  also  lower 
in  accordance  with  trend  iv  noted  above.    Altogether,  the  trans- 
ition observed  by  Keilin,  et  al.  ,  is  much  weaker  than  ours;  we 
find  from  their  report  that  their  wall  temperatures  were  no- 
where more  than  1  K  above  bulk  fluid  temperature  even  beyond 
the  transition.    It  appears  that  the  conditions  of  transition  at 
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Figure  16.    Heat  transfer  coefficient  as  a  function  of  quality  (comparison  of  this 
study's  data  with  that  of  reference  4). 
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L/ D  =  278   are  much  different  from   L/D  =  48  .    There  is  much 
more  vapor  present;  therefore  it  must  travel  at  higher  velocity 
and  the  heat  transfer  coefficient  after  transition  is  not  in  such 
contrast  with  that  before  transition,  which  is  the  case  in  our 
results , 

.  3  Comparison  of  Modes  of  Heat  Transfer 

In  figure  17  the  present  results  of  heat  flux   q  as  a  function 
of  wall  temperature  rise    AT   are  plotted  for   L/D  =  48  for 
comparison  with  other  modes  of  heat  transfer.    We  have  includ- 
ed calculated  values  of  the  nucleate  pool  boiling  correlations 
of  Kutateladze  mentioned  above  and  calculated  values  of   q  vs. 

AT   for  supercritical  helium  heat  transfer  at  3  atm  and  at 

4        5  6 

Reynolds  numbers  of  10  ,  10    and  10  .    These  were  computed 
from  the  correlation 


developed  in  reference  [  l]  which  represented  the  data  of  that 
reference    (L/D  =  20  and  40)    with  an  rms  percent  deviation  of 
8,  5%.    Supercritical  data  obtained  in  the  present  work  are 
represented  by  this  equation  to  within  10%  for  most  points  for 
which   L/D  >  20. 

We  conclude  from  this  comparison  that  for  heat  fluxes 

2 

between  0.01  and  0.3  W/cm    subcritical  pressures  are  prefer- 
able to  supercritical  pressures  for  equivalent  mass  velocities 
in  forced  convection.    Above  these  heat  fluxes,  for  the  data 
presented  here  supercritical  heat  transfer  is  superior  to  sub- 
critical,  which  is  now  in  a  dry-wall  regime  and  shows  very 
steep  temperature  rise  for  a  small  heat  flux  increment. 
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PERFORMANCE  OF  CENTRIFUGAL  PUMP  WHILE 
CIRCULATING  LIQUID  HELIUM 

During  the  performance  of  the  heat  transfer  experiments  there 
was  no  evidence  of  vapor  binding  or  instability  in  the  centrifugal 
pump  provided  the  loop  pressure  was  above  1.  1  atm.    The  pump 
appeared  to  perform  according  to  the  characteristics  shown  in 
figure  18.    These  performance  curves  were  determined  as  described 
in  reference  [2  ]  and  the  load  line  of  the  entire  flow  loop  was  calcu- 
lated from  standard  friction  pressure  drop  equations.    It  should  be 
noted  that  the  load  line  does  not  change  substantially  for  different 
densities  in  the  pump.    However,  knowledge  of  the  density  in  the 
pump  is  important  in  determining  the  mass  flow  rate. 

In  figure  19  the  flow  rate  measured  calorimetrically ,    m  , 

exp 

is  compared  with  that  derived  from  the  pump  characteristics  and 

the  calculated  load  line,    m  The  density  of  the  liquid  was 

caic 

assumed  to  be  that  at  the  heat  exchanger  temperature  and  the 

measured  loop  pressure.    The  abscissa  in  figure  19  is  the  subcool- 

ing  at  inlet  measured  in  watts.    For  subcritical  pressures  this  is 

given  by  m       (-^cat"  ~        where  H  is  the  enthalpy  of  satura- 

exp     k3 A  J.  jJ\.  i. 

ted  liquid  at  the  loop  pressure  and  H   is  the  actual  inlet  enthalpy. 
For  supercritical  pressures   ^g^-p        taken  to  be  the  enthalpy 
corresponding  to  loop  pressure  and  the  transposed  critical  tempera- 
ture.   It  is  seen  from  the  data  available  that  the  pump  performed  to 
within  ±  15%  of  its  predicted  performance  at  the  higher  pressures. 
The  poor  performance  at  1.  02  and  1.  097  atm  is  attributed  to  vapor 
formation  in  the  pump  due  to  heat  leak  down  the  pump  shaft.  This 
is  a  special  case  of  cavitation  which  is  entirely  dependent  on  the 
installation  of  the  pump.    We  estimate  a  heat  leak  between  0.  8  and 
4.  0  watts  in  the  present  installation  which  certainly  explains  the 
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Figure  18.    Pump  performance  characteristics  and  flow  loop  load  line. 
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Figure  19.     Ratio  of  measured  flow  rate  to  flow  rate  predicted  from  centrifugal 
pump  performance  characteristics  vs  subcooling  power. 
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poor  performance  at  1 .  02  and  1.  097  atm.    However,  at  higher  pres- 
sures the  pump  appeared  to  perform  in  spite  of  the  heat  leak  and 
this  is  not  entirely  understood, 

CONCLUSIONS 

i.  The  usual  transition  from  wetted-wall  to  dry-wall  heat 
transfer  is  observed  in  heat  transfer  to  subcritical  helium  I 
under  forced  convection, 

ii.  Below  the  transition,  forced  convection  has  little  effect  on 
heat  transfer  and  the  pool  boiling  correlation  of  Kutateladze 
for  the  heat  flux  vs,  temperature  rise  of  the  wall  is  suffi- 
cient. 

iii.  Above  the  transition,  the  heat  transfer  coefficient  falls  off 

drastically.    For  superconductivity  applications  this  will 

usually  have  to  be  avoided  except  at  mass  velocities  of  the 
2 

order  of  60  g/cm  s  or  higher. 

iv.  Transition  is  a  function  of  many  variables  (see  3.2).  More 
work  will  be  required  to  make  this  a  predictable  phenome- 
non, 

2 

V.     At  heat  fluxes  between  0.01  and  0.  3  W/cm  ,  subcritical 

heat  transfer  is  superior  to  supercritical  for  the  same 

2 

mass  velocity.    Above  about  0.  3  W/cm    the  reverse  is 
true, 

vi.    The  centrifugal  pump  used  to  circulate  the  liquid  performed 
according  to  the  characteristics  determined  earlier  pro- 
vided that  below  1.  1  atm  the  subcooling  power  was  not  less 
than  about  2  W, 
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